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1. Why this paper? 

Lately lattice QCD has enjoyed considerable success. In 2003 the HPQCD, MILC, and Fermi- 
lab Lattice Collaborations found that calculations based on Monte Carlo simulations agreed with 
experiment at the 1-3% level for a wide variety of physical quantities ||l|]. During the next two years 
the same techniques were used to predict the shape of the D — > Klv semileptonic form factor [§], 
the mass of the Be meson [^], and the decay constants of the D and mesons [Q]. Each of these 
results was subsequently confirmed by experimental measurements [^. Lattice QCD was also used 
to determine the strong coupling [^. These results have made a strong impression in the fields 
of particle physics and nuclear physics, and a long-standing claim of lattice gauge theorists is now 
generally accepted: the combination of numerical simulation and chiral perturbation theory [^] 
is a sound way to solve gauge theories, such as QCD. 

The key to this development was the incorporation of 2-1-1 flavors of sea quarks with the. fastest 
technique, "improved staggered fermions with the Asqtad action" [Q]. Unsurprisingly for compu- 
tational science, the fastest technique has some unresolved theoretical questions hanging over it. 
This state of affairs presents a great opportunity. It is simple common sense to reproduce the results 
of Refs. [|l], H, ^ with theoretically cleaner treatments of sea quarks. Such methods require 
more computation, so one can argue (persuasively) to devote more computer resources to lattice 
QCD than in the past. At the same time, there is a strong motivation (and obligation) to understand 
theoretically how, where, and why staggered fermions work so well. 

The crux of the unresolved issues is the fermion doubling problem. Staggered fermions cope 
with the problem partially, reducing the number of species from sixteen (in four space-time dimen- 
sions) to four. The empirically successful results cited above all use ensembles generated and made 
publicly available by the MILC Collaboration [|l^]. To reduce the number of fermion species (per 



flavor) from four to one, gauge fields in these ensembles have the partition function [11] 

1/4 



J q 



e ■^'5^"°% (1.1) 



where 5gauge is the action for the gluons, p^^^^ is the (improved) staggered discretization of the Dirac 
operator, and the subscript 4 on the determinant emphasizes the underlying number of species. The 
question, then, is whether [detn{p^^ + niq)]^ 1'^ is legitimate. 



This so-called "rooting procedure" is controversial. Critics seem to accept that Eq. (1.1) is 
valid in perturbation theory, where internal fermion loops are simply multiplied by i. Several 
interested and disinterested parties are investigating non-perturbative aspects, sometimes finding 
undesirable features though not, to my knowledge, fatal flaws. An apparent exception is the work 



of Michael Creutz [12, 13, 14, 15], including these Proceedings [16], which does claim that rooting 



fails. Although these papers have already been refuted []17p, the organizers of Lattice 2007 invited 
me to respond to Creutz's claims "on behalf of the staggered community." My charge is to "increase 
the confidence of the wider lattice community in the rooting procedure" {italics added). 

The organizers also asked me to comment on the relevance of lattice gauge theory to CERN's 
Large Hadron Collider (LHC). This request and the need to anticipate new critiques have shaped 
the organization of this talk. Sectionals an essay presenting the view that the LHC era will require 
fast lattice calculations and, thus, a broader understanding of validity of the rooting procedure. 
Sections |3| and ^ disentangle the main issues. Even without rooting, staggered fermions bring 
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in complications, which are reviewed in sect. |3|, focusing on a new quantum number "taste" that 
characterizes the remaining species doubling. Furthermore, even with full SU(4nj) taste symmetry, 
rooting has its peculiarities, which may be unfamiliar but should not be controversial; they are 
explained in sect. 0. These two strands are then brought together in sect. || to discuss rooted 



staggered fermions. It is then possible to refute Refs. [12, 13, n4, 15, 16] relatively tersely, in 



sect. ^ Section ^ covers some new developments, and sect. | offers a summary and perspective. 

2. The lattice community in the LHC era 

To represent the "staggered community" before the "wider lattice community" it helps to think 
about who belongs to each of them. The wider lattice community simply consists of those who 
attend conferences on lattice field theory and submit e-prints to the hep-lat section of the arXiv. The 
staggered community includes not only those who generate large ensembles of SU(3) gauge fields 
with rooted, staggered sea quarks, but also those (like me) who use them. The staggered community 
also reaches beyond hep-lat. Many physicists without much expertise in lattice gauge theory have a 



stake in the validity of the rooting procedure. For example, the Heavy Flavor Averaging Group [ |18[ ] 
and similar enterprises use lattice QCD results to help them understand whether non-Standard 
phenomena contribute to B-meson decays. It is therefore useful to take the staggered community 
to be everyone who thinks the validity of staggered fermions is worthy of study. 

The organizers' charge implies that there may be some in the wider lattice community who 
do not see themselves as part of the staggered community, even as defined here. That would be 
unfortunate. Despite the "controversial" nature of Eq. ([lTI), the arguments against it are difficult to 
pin down, because they are not documented in the scientific literature. Mike Creutz is to be com- 
mended for writing up his critique, making it possible to decide whether his qualms are correct. On 
the other hand, some members of the lattice community seem to prefer sniping from the sidelines. 
They may hope that increases in computer resources will make the controversy moot, by allowing 
other methods for sea quarks to catch up. The status of other methods is surveyed in Refs. [[T9|], and 
one may judge for oneself. My reading is that other methods have not caught up, partly because 
algorithms for the staggered sea are not standing still [pO|]. 

Future calculations are more saUent than past efforts, so it is worth looking ahead to the LHC 
era. The focus of the LHC's physics program will be on the terascale,^ where we expect to find the 
agents of electroweak symmetry breaking. They may be strongly coupled and, if so, the need for 
numerical lattice gauge theory will skyrocket. A lot will be at stake, and the fastest way to elucidate 
the physics will prevail. Even if physics at the terascale is weakly-coupled or, worse, boring, a need 
will remain for precision and accuracy in /k, fsV^^ ^^c- I^i addition, searches for non-Standard 
particles will profit from good calculations of moments of parton densities, especially the gluon 
density, so that signal and background cross sections can be calculated. 

In summary, whether for QCD or for other gauge theories, the fastest correct technique for 
simulating the fermion sea will remain a key tool during the LHC era. New researchers, young 
and old, may turn to lattice gauge theory, and they will expect that we established experts have 
understood whether rooted staggered fermions are valid, or not. For this reason I think the whole 



'Using the term "terascale" for the teraelectronvolt energy scales was inspired by terascale (i.e., teraflop/s and 
terabyte) computing (J.D. Lykken, private communication). 
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lattice community should engage in an open-minded scientific debate, and skeptics should submit 
their criticisms of rooted staggered fermions to the hep-lat arXiv. 



3. Staggered fermions without rooting 



Most of the complexity of rooted staggered quarks has nothing to do with the rooting proce- 
dure. Instead, it stems from the way spacetime and flavor-like symmetries of four (or 4n f) Dirac 
fermions emerge in the continuum limit. The aim of this section, therefore, is to give a brief review 
of the definitions and symmetries of staggered fermions and their connection with the continuum 
limit. Most of this is not rigorously proven, but it is fair to say that it is established, principally 
because numerical simulations shore up the theoretical picture. 

The simplest discretization of fermions replaces Dirac's covariant derivative with a nearest- 
neighbor interaction. The resulting "naive" action is 

5nai-ve = ^^^^^(x)/^ [U ^{x)T{x + fla) - U^x - fla)T{x - flu)] +moa4£r(x)T(x), (3.1) 

where each site x possesses Grassmann variables and T'^,, with / and a the color and spinor 
indices. This action is invariant under color SU(A'^t:), lattice translations, and hypercubic rotations. 
With Hf flavors there is a Uv(«/) x Ua(«/) chiral flavor symmetry, softly broken by the masses. 



The naive action also possesses a remarkable SU(4) "doubling symmetry" [|21|], with fifteen 
x-dependent (anti-Hermitian) generators B"^ (x) : 

B^IW = 7^I75(-l)"^ Bsix) = ijseix), B^{x)Bs{x), B^{x)By{x) < v), (3.2) 
where n=x/a, and e(x) = (_i)''i+"2+n3+"4 jj^g naive fermion field transforms as 



T(x) ^ e'»'«'Wr(x), T(x) ^ T(x)^'-'»'^''(-^) 



(3.3) 



The physical interpretation becomes clear in momentum space. Consider, for example, a doubling 
transformation generated by B^: T{p) coso)T{p) + smo)Y^Y5T{p + jxTi/a). In general, the 
doubling symmetries relate all 16 comers /aoi the Brillouin zone, up to a shuffling of the Dirac 



index. (See Eqs. ( 3.16 )-( 3.20 ), below, for a complete list of the 4- vectors n^.) 

The physical consequence of the doubling symmetry (in four dimensions) is that a single naive 
fermion field T corresponds to 16 species of fermion. The extra species are evident in vacuum 
polarization, leading to jSq = — |l6?iy in the running of the gauge coupling [^. The axial 



anomaly receives contributions from all 16 species, in the pattern [ ]21[ ] 



(0,0,0,0) 

1 



(;r,0,0,0) (;r,;r,0,0) 

4 + 6 



(;r,7r,;r,0) (K,n,n,n)' 
■ 4 + 1 



21 = 0, 



(3.4) 



where a typical tt^ hovers over the integer multiplicity of the species from that kind of comer. The 
total anomaly vanishes because with the naive action the flavor-singlet axial symmetry Ua(1) is 
exact. Because the anomaly is wrong, naive fermions do not seem to be what one wants for QCD. 



The doubling symmetries can be rendered x-independent via a change of variables [ 22 1 



r(x)=^l(x)x(x), T{x)=x{x)a-\x), a(x) = f,^i^il'f. 



(3.5) 
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Rewriting the naive action in the new fermion fields X'^,X'^, one finds 



\a^Y^X{x)r]^{x) [U^{x)X{x + jxa)-Ul{x-jxa)X{x-jxa)\ +moa'^£x(x)X(;c), (3.6) 



where the signs T]^ (x) have replaced the Dirac matrices via 

Qr\x)Y^Q.{x±fxa) = (-l)^p<f "p =: ^^{x). 



(3.7) 



Equation (3.6) can also be obtained by diagonalizing a maximal commuting subgroup of the dou- 
bling symmetry [23]. 

In Eq. (3.6) the transformed spinor index is sterile, so the number of degrees of freedom can 
be reduced four-fold, yielding 



?stag = \a^Y.^{x)r]^{x)[U^{x)x{x + (la) - Ul{x- fLa)x{x- (10)]+ niQa^ Y.^{x)x {x) , (3.8) 

X,^ X 



where %' is a fermion field without a spinor index. A Hamiltonian formalism with one-component 
fermions and sign factors instead of Dirac spinors and matrices was introduced by Susskind [^], 
extending work in l-i-l dimensions by Banks, Kogut, and Susskind [25]. Because of the ubiquitous 
factors (—1)"'' the Euclidean formulation with ^stag and is called "staggered fermions." 

The projection from the big to the little x' removes the SU(4) doubling symmetry, al- 
beit in a not-so-straightforward way. All other symmetries — except color and the vector flavor 
symmetries — ^become x-dependent [p4l p^. For example, consider translations 




J{x + (ia) 
T{x + (ia) 



X{x) ^ i;^{x)Y^X{x + (ia) 
X{x) ^ i;^{x)X{x + (la)Y^ 



where 



Q.-\x)Q.{x±(La)Yfi = (-1)^" 



(3.9) 



(3.10) 



Because of the Dirac matrix in ( |3.9[ ), we see that the projection from X'„ to x' does not commute 
with lattice translational symmetry. On the other hand, a certain combination of translations and 
B transformations called shifts does survive the projection: 



T{x) ^ -iB^B5T{x + (la) ^ + A«) 

t{x) ^ -it{x + (la)B^B5 x{x) i;,^{x)x{x + (Xa) 
Uv{x) h->-Uv{x + (la), Vv 



(3.11) 



and Xa transforms just like X- Similarly, rotations, spatial inversion, and charge conjugation be- 
come entangled with this residue of the doubling symmetries [^, 27, 28, 29, 30, pl| ]. 

Here we shall focus on issues related to the species content of staggered fermions. Acting on 
the fermion fields, the shifts anti-commute, S^Sy = SyS^ (v / /i), from which it follows that 
the residue of the doubling symmetry is a discrete Clifford group Because shifts translate the 
gauge field, a single r4 symmetry arises for any rif. Shift symmetry has two kinds of irreducible 
representations (irreps) — fermionic and bosonic — with representation matrices, respectively. 



(3.12) 
(3.13) 
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where the physical momentum ranges over G {—K/2a,K/2a]. The fermionic representation is 
4-dimensional, labeled by an index f = 1 , 2, 3, 4, and the are 4 x 4 matrices obeying the Clifford 
algebra, = 25^v- The sixteen 1-dimensional bosonic irreps are labeled by the corners A 

of the Brillouin zone and the pre-factors are signs, 5^ = e'^^' . The index t and the label A denote a 
quantum number that, nowadays, is called taste. 

A colored particle's taste, just like its momentum, is made complicated by SU(A^f ) color gauge 
symmetry. It is enough, however, to consider color singlets, such as the mixed-action bilinear, 
which arises, for example, in heavy-light physics p2[]. Let 

(3.14) 



where the matrix Q. reappears, now with the first index interpreted as a Dirac index, but the second 
index interpreted as taste index [F 



34]. Here the field *P represents a Wilson or Ginsparg-Wilson 



anti-fermion, so its spinor index has the conventional meaning. Under shift symmetry (as one can 

(r) 

easily verify), the meson field Hf ' transforms as 



: Hp (x) ^ [^^]n'H^,' > {x + fla) , ^^=y^ 



(r) 



(3.15) 



(r) 



so Eq. (3.15) gives an explicit realization of Eq. (3.12). Under rotations and spatial inversion 
transforms as a scalar, vector, or tensor, etc., dictated by F, up to a change of taste. Baryons with 
three staggered quarks also have fermionic taste [27,|35]. 

Staggered-staggered mesons have bosonic taste A. Rotations imply degeneracies in an almost 
obvious way,^ leading to multiplets 

/ : TT^ G {(0,0,0,0)}, (3.16) 

y : TT^ G {(71,0,0,0), (0,71,0,0), (0,0,7r,0);(0,0,0,7r)}, (3.17) 

r : 71^ G {(71,71,0,0), (71,0,71,0), (0,71,71, 0);(7r,0,0,7r), (0,71,0,71), (0,0,71,71)}, (3.18) 

A : 7r^ G {(0,71,71,71), (71,0,71,71), (7r,7r,0,7r);(7r,7r, 71,0)}, (3.19) 

P : 71^^ G {(71,71,71,71)}. (3.20) 

These multiplets are called the taste-singlet /, the vector taste V , the tensor taste T , the axial vector 
taste A, and the pseudoscalar taste P. Another hadron with bosonic taste is a heavy-light baryon 



consisting of two staggered light quarks and a Fermilab or NRQCD heavy quark [36]. 



The staggered action, Eq. (3.8), inherits softly broken but otherwise exact chiral symmetries. 



'•T"e(x) 



X{x), X{x)^x{x)e^ 



(3.21) 



where T"^ = —T". With Uf flavors, this is simply the Ua(«/) symmetry manifest in Eq. (0), 



made x-dependent via Eqs. (3.5). Thus, one has axial currents and pseudoscalar densities 



A];"(x) = \r]''{x)e{x) [x{x + va)Ul{x)T"x{x) - x{x)Uy{x)T'^x{x + va)\ 
P"p{x) = e{x)x{x)T-x{x), 



(3.22) 
(3.23) 



^Strictly speaking, the degeneracies are among irreps of the symmetry group of the transfer matrix, which entails cu- 
bic rotations, not hypercubic rotations [ ^jj ^o| ] . Thus, states with =Q,K need not be degenerate; hence the semicolons 
in Eqs. (5.17)-(3.19). In practice they turn out to be nearly degenerate [po[]. 
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satisfying PCAC relations [ ]37[ ] (for equal renormalized masses m)^ 



A^'p^x-va)] =2mP^{x)+2a-'^e{x)trT''. 



(3.24) 



It is tempting to say that the axial anomaly does not appear in Eq. ( |3.24D owing to Eq. ( |3.4j ), with 
each contribution at quarter strength. It is better, however, to note that it is not expected, because 



(3.25) 



the symmetry ( |3.21[ ) is exact. The axial anomaly is also undesired. Under shifts Aj," and Pp 
transform as 

, iA^p"{x) ^ {-lY^A^p"{x + fla) 
^ ■ [ P^{x) ^ i-\)"^'P^{x + ila) 

so one sees that these bilinears are not taste singlets, rather they transform under the (non-singlet) 
P irrep. This non-singlet character of Ap" and Pp provides a clue that the new quantum number 
taste plays a key role in the physics that emerges in the continuum limit. 



The change of variables ( |3.5[ ) followed by the projection from to x' reduces the number of 
species from 16 to 4. The action can be rewritten using physically suggestive fields 

¥at{y) = iY.^at{r)U{y,y + r)x{y + r), Yatiy) = ^'£xiy + r)^Ur)U{y + r,y), (3.26) 



where y labels hypercubes of size 2^ or, equivalently, sites on a coarser lattice of spacing b = 2a, 
r runs over the hypercube, and U{y,y + r) is parallel transport along some chosen path from y to 
y + r. As in Eq. ( |3.14 ), the indices a and t are interpreted as Dirac and taste indices. Suppressing 
the gauge field, the action is then rewritten [33, |4|] 



4 4 

?stag = WY.ll^' [V^' {y + f^b)-\ift{y-flb)]+ mob"^ Y^H^t (j) ¥t (y) 

y,^if=l y f=l 

4 

- i^^L L Vtiy)Y5[^5^ti]tt'[Vt'iy + iib) + Yt'iy-iib)-2YAy)], 

y\lit,t'=l 



(3.27) 



which looks like a lattice field theory of four Dirac fermions, with a Wilson-like term to alleviate 
doublers. For perturbative gauge fields this appearance continues to hold. For example, the beta- 



function starts with jSq = -jA'c — 3 4n/ [23], and the exact r4 part of the shift symmetry ensures that 



mass renormalization is taste-independent [26, 38]. 

The central conjecture of (unrooted) staggered fermions is that the picture of four (or 4nf) 
Dirac fermions holds non-perturbatively. I say "conjecture" to mathematicians because it is not 
proven rigorously (and I could say the same about all formulations of lattice fermions), and to 



physicists because the coupling to gauge fields is based on the one-component action, Eq. ( [3.8| ) 
Consequently the spacetime and flavor-like symmetries are entangled, so it is not especially trans- 
parent how S0(4) rotational symmetry and SUv(4?i/) x SUa(4?i/) x Uv(1) flavor-taste chiral sym- 
metry emerge in the continuum limit. 



To get an idea of the complications (see, e.g., Refs. [ |26[ , |30|] for details), let us consider some 
important cases. The hypercubic rotations (denoted SW4) are embedded 



SW4 C [S0(4) X S0(4)]diag C S0(4) x SUv(4) C S0(4) x SUv(4?i/) 



(3.28) 



^In the flavor singlet tiT" 7^ 0, but the last term vanishes in a average over a hypercube, so it is unimportant. 
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into the diagonal subgroup of the Euclidean group and a flavor-singlet subgroup of the SUv(4«/) 
group. The Clifford group abstracted from the shifts and the vector flavor symmetries are embedded 



r4 X SUv(?i/) C SUv(4?i/) 
in the vector flavor-taste symmetries. The axial flavor symmetries are embedded 

Ua(«/)cSUa(4«/) 



(3.29) 



(3.30) 



in the flavor-taste non-singlet axial symmetries. For sect, the embedding ( [3.30| ) along with 
Eqs. ( 3.22| )-( |l!25 ) are key. The exact chiral symmetries are not taste singlets. Therefore, they 
cannot have anything to do with anomalies, gauge-field topology, zero modes of the Dirac operator, 
or their consequences, such as 't Hooft vertices. 

If the emerging picture is correct, staggered fermions must possess a bilinear corresponding in 
the continuum limit to the flavor- and taste-singlet axial current. Analogously to the case of Wilson 
fermions, this is not a Noether current, but the explicit breaking is superficially of order a. The 
anomalous current extends over a whole hypercube [23, ^ ]: 



Af(x) 



32 



£ ri^'{x + c)T]i{x + b)T]2ix + b)T]3ix + b)T]4{x + b) [ 



b+c=d 



X{x + c + }la)Ul {x + c)U {x + c,x + b)x{x + b) 
X{x + b)U {x + b,x + c)U^ {x + c)x{x + c + jXa) 



(3.31) 



where d = (1 + 2 + 3+ 4)a, U{x + b,x + c) is the average of parallel transport over paths from 
x + b to x + c, and the sum is over all 4-dimensional diagonals of the hypercube at x. Under shift 
symmetry Aj is a taste singlet, as is the corresponding pseudoscalar density 



Piix) 



I 

16 



^ rii{x + b)ri2{x + b)ri3{x + b)ri4{x + b)xix + b)U{x + b,x + c)xix + c) 

b+c-d 



(3.32) 



A J and Pi satisfy a PC AC relation [ ]37| ] with the anomaly of four Dirac fermions [g3|]. Compared 



to Eq. ( |3.24| ) a complication is the need for renormalization [|37p, but this is analogous to the renor- 
malization of Ward-Takahashi identities of Wilson fermions. 

If staggered fermions do indeed provide a (4?iy )-species version of QCD, then general features 
of the hadron spectrum should provide clear numerical evidence. The chiral symmetries, and the 
standard line of reasoning for spontaneously broken symmetry, imply that the pseudoscalar meson 
masses should behave (with degenerate quark masses tUq) as 



= 2mqB + a^A^ + n^5ao5^j, 



(3.33) 



where a 

.1/4 



0, . . . ,njr — I labels flavor and t, G {P,A, T,V,I} labels taste. The dynamical quantities B, 
A^''', and n are of order A, the characteristic dynamical scale of the gauge fields. The exact (taste- 
nonsinglet) chiral symmetries of Eq. ( |3.21| ) imply Ap = 0. The "t]'-?!** sphtting" parameter /i^ can 
be generated only for the flavor- and taste-singlet meson. This pattern is shown (for = 2) in 
Fig. [l[ It is consistent with extensive numerical simulations in the quenched approximation, as 
well as with 2 and 2-1-1 flavors of (rooted, staggered) sea quarks [10, |^], including the scaling of 
the pseudoscalar splittings a^A^ [pO[]. 
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The approach to the continuum Umit can be clarified with two theoretical tools, the Symanzik 
effective field theory Q and chiral perturbation theory (;^PT) [41, 42, 43, 44]. In the Symanzik 
local effective Lagrangian (LE^), dimension-four interactions have S0(4) rotational symmetry 
and, by design, SUv(4?i/) x SUa(4?i/^) flavor-taste chiral symmetry, softly broken by mass terms. 
Dimension-six four-quark operators break these down to SW4 and r4 x SUv(«/) x Ua(«/), as 
indicated by the embeddings ( 3.28 )-( 3.30 ) 42, 43]. In particular, some of these operators 
are invariant only under the r4 taste symmetry, not a full SU(4). They are variously called taste- 
breaking or taste-exchange interactions. The LE^ analysis clarifies why smearing strategies re- 
duce the strength of taste-symmetry breaking [45]. If successful, smeared actions should yield 



smaller splittings A^. They do []46[]. They also improve the scaling of the p -meson mass and the 
static potential [0]. 

The splittings and other effects of broken taste symmetry are, of course, a complication. They 
lead to multi-parameter fitting procedures, for example of the pseudoscalar decay constants [p^]. 
Another example is the n-n threshold, which splits into five n^-n^ thresholds, ^ G {P,A,T,V,I}. 
In last year's plenary talk on this subject [p9|], the need for complicating fitting was one reason to 
deem lattice QCD with rooted staggered quarks "ugly." Many of the complications come not from 
rooting but from the intricate symmetry structure and the desire for statistically sound fits. Fitting is 
nicely illustrated by /;;: and f^. The staggered data start so close to the experimental result that just 
about any chiral extrapolation would agree with experiment, but statistically good fits are obtained 
only when staggered ^PT is used [|^. Thus, these fits provide further evidence that staggered 
(valence) fermions do indeed simulate 4n f flavors of quark. 

To study how these properties of staggered fermions relate to gauge-field topology, one must 
consider the eigenvalues of the staggered Dirac operator. The exact (non-singlet) chiral symmetries 
imply that eigenvalues come in complex conjugate pairs, it/A,- + mq, with orthonormal eigenvectors 
fi{x) and e{x)fi{x). In a gauge field with a non-zero topological charge Q, the continuum Dirac 



m 




1= 1 



7 = 



Figure 1: Pseudoscalar spectrum for 2 flavors of staggered fermion, so 8 species in all. The isovector 
multiplets (with = +1, 0, —1) each consist of sixteen states, split by lattice artifacts of order A^a^ into 
submultiplets with 1, 4, 6, 4, and 1 states (for, respectively, irreps of taste I, V, T, A, and P). The isosinglet 
multiplet is similar, except that the taste-singlet I splits from the others by continuum-QCD effects (as usual). 
This 77'-like state suffers from noisy correlators, and numerical data are consistent with this picture without 
being definitive 
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operator has n± zero modes of chirality ±1, satisfying the index theorem — n^Q, where rig 
is the total number of species. With lattice staggered fermions, the zero modes are no longer exact. 
Nevertheless, if staggered fermions do indeed yield four quarks, then several features must emerge 



dynamically [pl[]. First, low-lying eigenvalues must cluster in quartets. Second, some (pairs of) 
eigenvalues should be near-zero modes with A,- ~ A^a. Third, the chirality of the corresponding 
eigenvectors should be close to ±1, in a way that satisfies the index theorem. Finally, in the finite- 
volume e regime, the distributions of eigenvalues should agree with expectations from random 



matrix theory [ 52 1 



Early numerical work was inconclusive [53], but with improved staggered actions, quartets 



of eigenvalues clearly emerge in numerical simulations [54]. In particular, if the gauge field has 
topological charge Q, then the eigenvalue spectrum contains quartets of the same (taste-singlet. 



Eq. ( |3.32D ) chirality, such that the index theorem is satisfied []54|]. Even with smearing, early 
comparisons of random matrix theory suggested that staggered fermions are topology blind [F5]. 



Equation ( |3.33| ) implies two regimes of interest [^g, |57|], however. 



£ regime : A » L ' » , Vi^ , 
e' regime : m^_^p » L^^ » mp. 



(3.34) 
(3.35) 



In the £ regime all "pions" are pseudo-Goldstone bosons, whereas in the e' regime most of them are 
very massive particles. In the £ regime staggered fermions should and do behave like a 4-species 



theory [57]; in the e' regime they should not [ pq ] and do not [55]. 

In summary, naive fermions appear to be problematic, because the anomaly is not generated. 
Analysis of the unexpected doubling symmetries rescues the formulation, however, yielding stag- 
gered fermions. Spacetime and flavor-like symmetries become entangled, making the interpretation 
as Dirac fermions less transparent. But now not only are the exact chiral symmetries non-singlets 
under a flavor-like quantum number called taste, but also an anomalous taste-singlet axial current 
can be found. On this basis, the emergence of 4n f Dirac fermions is theoretically plausible. Owing 
to a wealth of results examining the nonperturbative content and structure, it is also fair to say that 
the validity of (unrooted) staggered fermions has been established numerically. 



4. Rooting with full SU(4) taste symmetry 

The previous section reviewed some aspects of staggered fermions, focusing on how four 
species emerge in the continuum limit. Before discussing whether the fourth-root procedure can 
reduce these four species to one, I would like to present a Gedanken algorithm. The aim is to sepa- 
rate some potentially confusing aspects of rooting, free of the complications of staggered fermions' 
taste-exchange interactions. 

Suppose that an algorithm designer with a good imagination (and a wicked sense of humor) 
found a way to speed up "your favorite lattice fermions" by substituting 

det 1 {p+ m) = {det4 [{p+m)(^U]Yl\ (4. 1) 

thereby introducing four "tastes." Here deti is a determinant for rif flavors, with I taste per flavor; 
det4 is for n f flavors, but 4 tastes per flavor. If the determinant is real and positive, this step does not 
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change the Monte Carlo weight at all, because the right-hand side is just a clever trick to calculate 
the left-hand side. Of course, the trick fails if the left-hand side is not real {e.g., for nonzero 
chemical potential [58]), or if it can be negative {e.g., m < 0). 

To probe the dynamics of this system, first introduce sources {J",J^) for scalar and pseu- 
doscalar meson operators {—i\lfT"\if, —iYT"j5W) the original theory: 



Z[r,Jl] = J ^^{deti(0+m + 7 + /575)}^'- 

= y" W {det4[(^+m + 7 + 7575)® I4]} 



Nr/4 



(4.2) 
(4.3) 



where is the gauge-field measure (including e "^sange)^ and the number of "replicas" A^,- will be 
useful below. The T" are flavor matrices (with r° = il / syiuf), and 7(5) : 



Spontaneous symmetry breaking is revealed by looking at the Legendre effective action 



J[i)=J{i){o,it) 

where 7(5) (a, tt) are defined implicitly, so that they create prescribed fields (a", tt") 



(j"{x) = —- 



ZdJ"{x) 



'(5)-J(5) 



^ ^ dZ 



ZdJIix) 



(4.4) 



(4.5) 



and (|4.3[) as the invariance of the determinants under a change of basis. One finds [ 59 1 



The mass matrices for a" and 71" are obtained from second derivatives of P. In the case at hand, 
they are constrained by symmetry — invariance of the fermion action is expressed in Eqs. (51 



(4.6) 
(4.7) 



+ 



jCdb^d 



0, 



oc mass and anomaly terms, 



where the vacuum fields (acTTo) are those in the absence of sources 7(5).^ If one assumes that all 
vacuum fields vanish except the flavor-singlet scalar (7q , then one obtains a constraint 



oc mass and anomaly terms, 



(4.8) 



and no constraint on d^T/da^da^ . Equation ( p^ ) is the basis for formulae such as Eq. ( 3.33 ). Of 
course, the dynamical assumption that Oq has not been proven mathematically but has been 
established numerically with various types of lattice fermions. 

"^)— SU(4«/) X SU(4?i/) instead of 



Equation ( |4.3[ ) apparently has more symmetry than Eq. 
S\J{nf) X SU(?i/). I will call the extended symmetry a phantom symmetry, because it is a figment 
of an algorithm designer's imagination. One can study the dynamical consequences of the phantom 
symmetry by promoting the sources to taste-nonsinglets [ p7| ] : 

(4.9) 

^Instead of fixing 7(5) in Eq. ( 1-.5 ) to get prescribed fields, simply set 7(5) = on the right-hand side to get (ctq, Sg) 



Z[7^,7^]= j {det4[(^+m)®l4 + 7 + 75 75]} 
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where the flavor-taste generator index A now runs from to {4nf)^ — 1. One finds the same pattern 
of spontaneous symmetry breaking as in Eqs. (p^-(^.7), leading to mass relations like (O) also 



for the taste nonsinglet phantom pseudoscalars. In all, there are {4nf) — 1 pseudo-Goldstone 
bosons, instead of just nj- — 1, and one taste- and flavor-singlet T]'-like meson. 

But what should one make of the phantoms? The equivalence of the two formulations, ( p!2[ ) 
and (pT3|), shows that phantoms do not exist. The symmetry analysis suggests they do. (A numerical 
simulation of your favorite rooted fermions could look at taste-nonsinglet correlators to decide and 
would probably find them.) The conundrum is easily resolved if the phantoms violate unitarity and 



cancel each other out in correlators that are oblivious to taste [ 17 ] . The lack of unitarity, particularly 
when constrained by symmetry, is not a concern and, in fact, is familiar in gauge theories [|60|]. 

The Gedanken algorithm helps establish a foundation for some questions concerning rooted 
staggered fermions. First, we see what kinds of correlators are physical. For bosonic taste, these 
are taste singlets and anything related to taste singlets by symmetry.^ For fermionic taste, one 
needs single-taste correlators [p9|]. It is possible to construct unphysical correlators, which could 
lead to confusion or mistakes. Another use of Eqs. (4^) and ( p~3| ) is to set a criterion for proving 



staggered fermions incorrect. Any argument that would also kill the Gedanken algorithm simply 
does not pertain to the issues at hand. An example is the proof that {p+my/'^ is non-local [jsl]], 
which does not prove that [det4(^+m)]'/'* cannot be obtained from a path integral with a local 
action []62|]. Attempts to disprove staggered fermions must focus on features not shared with the 
Gedanken algorithm, namely the violations of perfect SU(4) taste symmetry. Finally, note that the 



details of the phantom sector depend on A'^r- For Ny G the rooted formulation of Eq. (1.3) is 



identical to a manifestly sound quantum field theory. On the other hand, Eqs. (4-.2) and (4.3) show 
that irrational choices, Uke Nr = 71, are simply irrational. 

5. Rooting with staggered fermions 

We are now ready to discuss rooted staggered fermions, which is the Ansatz that quarks can 
be simulated with Eq. ( |1.1[ ), using any of several available algorithms [ ]63| , An initial set of 
arguments in favor are as follows. First, Eq. ( |1.1[ ) resembles Eq. (p^), especially when one looks 



at the ^jjjjg defined via Eq. ( |3.27[ ); the four-taste determinant has SU(4) taste- violating parts that are 



suppressed superficially by powers of the lattice spacing. Second, perturbative renormalization [ ]26| , 
38] and the nonperturbative features discussed at the end of sect. |3| support the picture of unrooted 
staggered fermions as a QCD-like theory, in which violations of SU(4) taste symmetry vanish in 
the continuum limit. Finally, having thus achieved a 4 x 4 block structure, the fourth root is as 
straightforward as in sect. ^ 

Implicit in this line of reasoning is the assumption that the continuum limit and rooting com- 
mute. We shall come back to this below. First, however, I would like to consider objections that 
hold even if these arguments are all correct. In particular, I stipulate for now that SU(4mj) flavor- 
taste symmetry emerges in the continuum limit, citing the evidence presented at the end of sect. ||. 

In the context of rooting most of the SU(4?i y) symmetry is a phantom symmetry. As such one 
should ask whether a phantom sector of particles is seen in numerical simulation. (They are, in the 
scalar propagator, see below.) If so, they violate unitarity. Moreover, at non-zero lattice spacing the 

^For example, one may use the taste-P irrep for pseudoscalar mesons, the taste-V irrep for vector mesons, etc. 
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symmetry is smaller than SU(4?iy), as exhibited in (3.29) and (3.30). Consequently, the phantoms 
are not degenerate, so they will not cancel identically in physical correlators. Fortunately, sect. § 
shows in general (if sketchy) terms that the rooted dynamics provide a safe house for phantoms, 
when (as stipulated above) they become degenerate in the continuum limit. Then the phantom 
sector is benign, even if it is not unitary, because the cancellation becomes perfect. 

The ramifications of non-degenerate phantoms are still ugly, in some cases perhaps even scary. 
A useful and illuminating tool is ;^PT adapted for the case at hand: rooted, staggered chiral pertur- 
bation theory (RS;^PT) ||6^]. This extension of partially quenched ;^PT encodes the A^,- dependence, 
so it is sensitive to the dynamics of the sea. As mentioned in sect. ^, the fit to pseudoscalar meson 
masses and decay constants provides clear evidence that this description works. To isolate the sen- 
sitivity of the numerical data to the sea quarks, MILC now tries fits in which the number of replicas 
floats. At this conference Bernard reported [65] (in his notation, Nr = 4nr) 

lN,. = nr = 0.2S{2){3), (5.1) 

in striking agreement with the expectation 1/4. Further evidence comes from heavy-light decay 
constants. Continuum partially quenched ;^PT blatantly fails to describe the numerical data [Q], 
while RS;tPT yields a statistically sound fit [Q]. The topological susceptibility is also better- 
described with RS;^PT than with continuum ;^PT []67|]. 

The pseudoscalar mesons just mentioned are low-lying states, but higher mass particles must 
be considered too. Correlators, for example those of glueballs, contain not only the single hadron of 
interest but also two-phantom-particle states. It can become hard to determine the single-particle 
mass, but these correlators are still interesting, because the two-particle contribution probes the 
sea []68|]. Correlators of taste-singlet scalar mesons contain Tl^-rj^ pairs, = P,A,T,V, and 
and T]^ consist of a valence and a sea quark. These Tl^-rj^ bubbles can be described with ;tPT, 
analogously to the quenched scalar correlator [|^]. The bubbles' weights depend on Nr- For A^,- < 4 
some weights are negative, as one expects when unitarity is violated. Once all thresholds become 
degenerate, the bubbles always add up to the correct total. Comparisons of numerical data with the 
RS;tPT formulae confirm this picture [70|, supporting the correctness of the rooted sea. A similar 
story holds in the case of a mixed action with a rooted sea and Ginsparg-Wilson valence quarks; 
the numerical data are again well described by mixed-action ;^PT Jti]]. 

In summary, the violations of unitarity are a serious objection and can be a practical obsta- 
cle. The numerical data suggest, however, that these effects, including their influence on the error 
budget, can be managed with RS;^PT. Furthermore, the analysis of sect. ^ reveals a safe house for 
phantoms, not only for pseudoscalar mesons, but also for other hadrons. 

A distinct, but related, issue is the locality of rooted staggered fermions. Critics of the rooting 
procedure have long suspected non-local behavior [72], but to my knowledge the concern is that a 
Lagrangian of the form YiP^ m)^/^!// is non-local [61]. Because this line of attack would kill the 
manifestly correct Gedanken algorithm, this concern is a red herring. There is, however, another 
kind of non-local behavior. Suppose one writes 



det4 (^stag + m)= det 1 (^u(4) + 



4 



(5.2) 



where the first factor on the right-hand side is manifestly SU(4) symmetric. The remaining factor T 
cannot represent a set of local interactions for the gauge field, because the left-hand side generates 
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a spectrum with taste splittings, and the SU(4)-symmetric determinant does not [73]. It is certainly 
accurate to call this a non-locality, but I am not sure how illuminating it is. Before taking the fourth 
root, the non-locality is particle-like. After taking the fourth root on both sides of Eq. (^^, it 
strikes me as plausible that the non-locality of T^^^ does not ruin cluster decomposition, it may 
become local in the continuum limit [74|, and it certainly is not the non-locality of a propagator 
(p+m)^^^^. Thus, it is not the non-locality that critics seem to have in mind. 

To be confident that the continuum limit is correct, it would be helpful if one could establish 
a Symanzik LEJf , not with an Ansatz based on symmetries, but through a derivation from the 
underlying lattice field theory. The Symanzik structure is evident at the tree level, and it seems 
to be accepted without controversy at all orders in perturbation theory. (It is not, however, proven 
at the level needed to prove Reisz's theorems [75]. An important first step has been to extend the 
power-counting theorem to staggered fermions [76].) In a forthcoming paper, Bernard, Golterman, 
and Shamir (BGS) show a new way to do so. For a short version with more details than I can give 
here, see Golterman's talk at this conference [77]. 



The new work of BGS is based on Shamir's block-spin renormalization group [pSQ, which was 
covered thoroughly in Sharpe's review [ p9| ] last year. After n blocking steps, Shamir arrives at a 
blocked staggered operator 

{p+m)„(E)U + aAn (5.3) 

with a SU(4)-symmetric part {p+m)n 14 and a taste-breaking defect A„. The blocked quantities 
live on a lattice with spacing 2"a, and in the continuum limit n — > 00, a ^ 0, cif = 2"a fixed. 
A power-law divergence in A„ is expected to be a^^, implying that the combination aA,, vanishes 
in the continuum limit. 



The blocked determinant resembles Eq. ( [4.9| ), with sources set to provide aA„, and the 
complicated measure for the gauge field in Ref. [78]. BGS develop an expansion in the underlying 
spacing a with the aid of some bookkeeping devices. The determinant is generalized to be 

s 1 . A/ {det4 [(£)+ m)„ (g) I4 + faA„l I"-- 

where A^^ need not equal A^^ = An^, and t need not equal 1. Sources for valence fermions are 

expjf] [(|>+m)„(g)l4 + vaA„]"\|, (5.5) 

where, again, for bookkeeping v need not equal 1. Post-analysis, one may set A^^. = A^,-, t = v = 1. 

Expressions (5.4) and (5.5) can be expanded in t and v, justified by the small quantity aA„. The 
double-expansion can then be reverse-engineered to reproduce the Lee-Sharpe LE^ for staggered 
fermions [pT]]. Each bookkeeping parameter brings an advantage. Possible nonperturbative non- 
polynomial dependence of the LE^ on A^^ is kept, because other factors are expanded in t and v. 
The expansion in t ensures that, to any order in a, the dependence on n,- is polynomial, so analyti- 
cal continuation from an integer to 1/4 is allowed. The expansion in v shows that valence quarks 
control the symmetries and the field content of the LE^, which do not depend on rir. Pending 
confirmation of the assumptions built into the blocking [78, |^], these are a strong results. 

This analysis underscores the importance of checking numerically that A„ scales in such a way 
as to justify the expansions. There are several pieces of evidence to suggest it does. Exhibit A 
is all experience with anomalous dimensions in QCD. Exhibit B is all experience with the pseu- 
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doscalar spectrum as taste-breaking of valence fermions is reduced [46, |T^; heuristically this is 
like reducing v without reducing a. Lastly, Exhibit C is a pilot investigation of A„ itself [79|. 

A lack of scaling of aA„ would kill rooted staggered fermions, because then we would be left 
with a non-unitary, non-local theory, even in the continuum limit. It is therefore worth stating how 



the evidence given above could be misleading. The pilot investigation [ J79[ ] is not yet definitive. The 
pseudoscalar splittings probe only four of the plethora of taste-exchange effects generated by A„. 
Finally, it could be that rooting and the continuum limit do not commute, in a way so profound that 
the rooted determinant itself generates A„'s anomalously large anomalous dimension [p9|]. Such 
"self-inconsistency" is possible, but highly implausible. 



6. Explicit refutation of Refs. JJ, 0, |6|1 



With the preceding sections' outline of rooted staggered quarks as a basis for discussion, I can 
now address Mike Creutz's specific criticisms, which below are summarized in slanted font. Many 
of my counterarguments are the same as in the original refutation |[r 



6.1 Order of limits 

Rooted staggered fermions require a "peculiar" order of limits, a — > with m fixed, followed by 



m^O / |i2| , I J 4 |i4 /. The required order is especially peculiar with one flavor 

For two or more flavors, this criticism is misguided. Computers have a finite memory, so 
one takes the continuum limit, a — > 0, at fixed spatial volume, Lr'. But there is no spontaneous 
symmetry breaking in a finite volume. To select a finite-volume vacuum close to the infinite- 
volume spontaneously-broken vacuum, explicit symmetry breaking is needed. For this general 
reason, one should keep m / while taking the continuum limit, then take L ^ oo, and last m — > 0. 

Also, the continuum limit must be carefully specified whenever exact symmetries imply pseudo- 
Goldstone bosons [pO|]. Suppose there are two kinds of particles with masses m\a^ = kL, m^a^ = 
I?, and consider an unconventional family of continuum limits, k = niqa^^'', £ = Aa'^^. Then 



m 



I = niqA, = A^a^^P, so if > then tt's correlation length diverges while a's does not. This 
is a different universality class than the standard one, p = 0. Staggered fermions (without rooting) 
are subtler still, because the would-be Ti-like particles, except the one with taste P, have a mass 
rn^^a^ = + m| = m^A + A^a^^^^. For p > 1/2, the continuum limit strands these at the 
cutoff, reminiscent of the e' regime. 

With one flavor there is no spontaneous symmetry breaking — the lowest-lying pseudoscalar 
is an T]'-like meson. In this case, the above considerations no longer apply. But every numerical 
lattice QCD calculation has unphysically large quark masses, to allow the algorithms to run faster, 
so I do not consider the order of limits to be a serious criticism even with one flavor. 

6.2 Mutilated quark-mass dependence 

Rooted staggered fermions yield the same system for —m as for m; therefore, the small-mass 
behavior is a function of np-, which we know is wrong in specific cases [12, 13, 14]. 



Odd powers of the quark mass stem from zero modes of p. As is obvious from Eq. (11), the 
rooting procedure turns m into (m^)^/^ = |m|. This defect of the algorithm has nothing to do with 



staggered fermions [17]. 
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It is instructive to examine how odd powers arise when one has near-zero modes A, ~ A^a, as 



with staggered fermions. Two pairs with eigenvalue ib/c,A^a, / = 1,2 arise [54], so 



1/4 



det4(0,„g + m) ■ oc [(c?AV + m2) {4A^a^ + m^)] ^'^ w |m 



(6.1) 



when the limits are taken in the correct order. If one takes the limits in the wrong order, then the 
determinant is indeed a function of n^. It is correct to consider this the wrong theory, but it is 
incorrect to assert that the numerical work takes the limit in the wrong order. The applicability of 
Eq. ( |6J] ) has been demonstrated clearly in the Schwinger model [81|. It is also possible to follow 
the m dependence of observables down until m ~ A^a^ <C A^a using RS;|^PT [f2|], well beyond the 



regime where the approximation in Eq. (6.1) applies 



6.3 Cancellation of non-unitary contributions 



The cancellations among taste multiplets seem contrived i[JJ, |J3[ |J4| ]. 

The Gedanken algorithm of sect. ^ shows that phantoms cancel each other automatically in the 
absence of taste-exchange processes ^VT\\. Thus, if the Shamir defect scales so that one can treat 
it as a perturbation, it is not reasonable to call the cancellation contrived. The scalar propagator 
provides a good numerical test, for two reasons. First, the weights of the non-unitary contributions 
depend (in ^PT) on Nr in an illuminating way, which could be checked numerically with (some) 
new simulations with Nr = 1,2,3,4. Second, the cancellation of the phantom modes depends on a^, 



which is being monitored as part of ongoing simulations [|70|]. 



6.4 Rank of chiral symmetry 



The rank (i.e., number of diagonal generators) of the chiral symmetry is wrong / |i4 |J4| ]. 

With four tastes the rank is indeed larger. But the Gedanken algorithm shows that this is not 
a feature of staggered fermions per se, but of the rooting procedure itself. It provides a clear ex- 
planation of why the physical sector is the taste-singlet sector. The other particles are phantoms. 
In particular the extra neutral pseudoscalars cancel each other out, perfectly in the Gedanken algo- 
rithm and (on the basis of numerical results) to order a^ with rooted staggered fermions. 

6.5 Anomalies 

The conventional axial anomaly cancels [12, 

This claim is simply wrong. The taste-P axial current suffers no anomaly, as desired, because 
it is off-diagonal in taste. The conventional staggered anomaly appears in the taste-singlet PCAC 



relation, as explained in the discussion of Eqs. ( [3.3 1[ ) and ( [3.32[ ). It comes with the right strength, 
because rooting multiplies the 4-species anomaly with the factor ^ appropriate for one species. 

6.6 Topology 

Rooted staggered quarks implement topology incorrectly, because quartets must break up at the 
boundary of topological sectors J^. Rooting averages over positive and negative chirality 

modes, so the index theorem cannot be satisfied /ll^ 



Studies of the eigenvalues reveal a quartet structure, including quartets of near-zero modes [ 54 ] 
As the gauge field transits from one topological sector to another, this structure must indeed be dis- 



16 



Staggered fermions: why not! 



Andreas S. Rronfeld 



rupted. But gauge fields near the boundary of topological sectors should have a large gluon action 
and are, hence, suppressed. If not, it is a drawback of the gluon action and has nothing to do with 
staggered fermions. It would be interesting to monitor this with smooth gauge fields. It would 
bode well for rooting if the quartets rearrange themselves rapidly, similarly to the way a lattice 
approximant to Q does. 



The assertion about positive and negative chirality is false. Ref. [57] shows that near-zero 



modes appear in quartets of the same taste-singlet chirality, Eq. ( |3.32D . The confusion may stem, 
as in sect. 6.5, from contemplating taste-P currents and densities. But taste-P chirality vanishes for 
all modes, including near-zero modes: Y,xfi{x)^ix)fi{^) = 0> because if fi{x) is an eigenvector of 
eigenvalue iXi + m, then e{x)fi{x) is the eigenvector of —iXi + m, and they are orthogonal. Hence, 
many statements about chirality in Refs. [12, [l^, 14, 15, 16] are simply ill-conceived. 



6.7 't Hooft vertices 

't Hooft vertices generate contributions to correlation functions that diverge in the chiral limit as a 
power of m^^ [15, It 



Only a full analysis of 't Hooft vertices can refute this assertion, because no explicit equations 



are given in Refs. [lO, ITH]. Here I shall examine the 77 '-like meson with one quark flavor, which 



should \ \15[ |16| ] expose the problem. The taste-singlet propagator consists of two terms 

C{x,y) = tr[G{x,y)U5G{y,x)U5], 
D{x,y) = tr[G(x,x)[/5]tr[G(3;,3;)[/5], 



(6.2) 
(6.3) 



where the trace is over color, and U5 abbreviates the sign factors, link matrices, and translations 
indicated in Eq. ( |3.32 ). The quark propagator G{x,y) = {x{^)x{y))x,x- ^^^^^ ^'^ combine 
C{x,y) and D{x,y) is 

{7l'{x)7]'{y)) = {-lC{x,y) + i-,D{x,y))^ , (6.4) 
where the sign arises from Fermi statistics. The weights ^ and ^ are crucial and follow immedi- 



ately from Eqs. ( |4.2| ) and ([1.3|). 

't Hooft vertices arise from (near) zero modes. At first glance, such gauge fields are suppressed 
by small eigenvalues from the determinant. For rooted staggered fermions the determinant factor is 



that given in Eq. (6.1). In a correlation function, however, these factors can be cancelled by small 
eigenvalues in the denominator, coming from fermion propagators. 



Let us consider the 2=1 sector for simplicity. Numerical simulation [ 54 1 tells us there are four 
near-zero modes, which we shall label ±1, ±2. Inserting an eigenvector-eigenvalue representation 
of the fermion propagators, the disconnected contribution is 



isolating the most singular parts as m 



7=±1,±2 



m 



f]{y)U5fj{y) + - 



(6.5) 



(with A <ti m). Similarly, the connected correlator is 



C(x,j)=It 



1 



1 



iXj + m 



fn^)U5fj{x)^^—-f]{y)U,fi{y). 



(6.6) 



Within a quartet, it is plausible to assume fj {x)Uifj{x 
should have a different taste. Then, 



iXj + m 

0(a), / / j, because each eigenvector 
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Ci^^y) = L ( -T^) fIi^)U5Mx)fJ{y)UsMy) + . . . ~ ^. (6.7) 



(=±1.±2 



When combining the two pieces according to Eq. ( |6.4| ), no singular behavior appears as m ^ (af- 
ter a —>■ 0). Thus, modulo one easy-to-check assumption, Creutz's claims about the 't Hooft vertex 
do not hold up in this, the simplest, example. Similarly, the 't Hooft vertex of the physical eight- 
fermion operator discussed in Ref. [Ill] is not singular, when Fermi statistics and orthogonality in 



taste is taken into account [83]. 



6.8 Summary 



The criticisms of Refs. [[IJ, [13|, [14j, |15|, |16p are based on three sources of confusion. Difficulties 
with mass dependence, the rank of the flavor-taste symmetry, and the way that phantoms cancel 
are explained by the Gedanken algorithm of sect. |^ This approach also illuminates the blunder 
of confusing the exact taste-nonsinglet chirality with the conventional taste-singlet chirality, which 
connects correctly to topology and 't Hooft vertices. Finally, it is easy to draw incorrect conclusions 
by choosing the incorrect order of limits (obdurately forcing m ^ before a ^ 0). The correct 



order is not "absurd" [12, 14], but necessary whether or not staggered fermions are employed. 



7. New developments 

There are two noteworthy methodological developments in improved actions for staggered 
fermions. One is a new discretization reducing taste-exchange interactions, while maintaining 
O(a^) improvement. The other is the completion of the full O(ava^) corrections to the gluon action. 
Either or both could be incorporated into future simulations of sea quarks. 

7.1 Highly improved staggered quarks 

Although the Asqtad action has much smaller taste-changing effects than the original staggered 
action, it would be better to reduce them further. The Asqtad action is obtained from the standard 



staggered action, Eq. ( [3.8| ), in two steps, as follows. First smeared links are constructed, 

where Ap is a covariant nearest-neighbor second derivative, and the product is symmetrized over 
all possible orderings of the directions orthogonal to pt. It yields bent staples of length 3, 5, and 7; 



substituting V for U in Eq. ( [3.8| ) yields the FAT7 action. The smearing is designed for, and is 
successful at, reducing the size of taste-exchange interactions. It does not achieve Symanzik im- 
provement, however. This is achieved by adding two improvement terms to the FAT7 action, the 
Naik [|^] term and the Lepage term [||], to obtain the a-squared (Asq) action and, with tadpole- 
improved couplings, the Asqtad action. 



A new action [ ]85[ ] extends the philosophy behind the FAT7 and Asqtad actions, introducing 
= ^ii'^ -^^U^, where the operator ^ brings the smeared V-links back into U(3), thereby 
reducing ultraviolet fluctuations, including taste-changing interactions [p^]. Using W in Eq. (Ob 



and generalizing the Naik and Lepage terms yields the highly-improved staggered quark (HISQ) 



action [85]. A numerical implementation [B7| is only around two times slower than Asqtad. 
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The first results from the HISQ action shed light on the issues of sect. Bl Reference [85] 



calculates in perturbation theory the size of the four-quark interactions in the Symanzik LE^ for 



both Asqtad and HISQ, finding the latter's to be an order of magnitude smaller. Reference [ |85[ ] 
also computes the a^A^ splittings of Eq. ( 3.33 ), finding them to be half as large with HISQ. In the 



language of and (^^), the HISQ action reduces the size of the (valence) Shamir defect, vaA„, 
by methodological means, and yields a valence spectrum closer to that of the continuum hmit. 

7.2 nf dependence of the gauge action 



The MILC ensembles [10] use the order-c? improved gauge action [88]. The pure-gauge 



one-loop matching of this action has been available for a long time [|89|]. Using these improved 
couplings, as MILC does, removes errors formally of order agNcCp' but not astifo^. Earlier this year, 
Hao et al. completed the fermion loop calculation with Asqtad fermions [pO|]. They find that the 
fermion loop has the opposite sign from the gluon-i-ghost loops. For = 3 they change the sign of 
the radiative correction to the coupling of the rectangle and reduce greatly the radiative correction 
to the twisted parallelogram. They recommend putting these results into future simulations and 



note that the size of the effects is what is needed to explain observed scaling violations [91]. The 



MILC Collaboration plans to use this result in future simulations [92] 



8. Conclusions 

Staggered fermions are fast, but not easy. After trivializing the spinor index of the naive 
fermion field, the projection to one component entangles flavor symmetries with spacetime sym- 
metries, and the remaining species doubling is reflected in a new flavor-like quantum number, 
taste. The exact flavor and taste symmetries are expected to enlarge to (softly broken) SUv(4«/) x 
SUa(4?i/), a mechanism that has been established by numerical simulations. 

When the fourth root is taken, staggered fermions have both too little symmetry (discretization 
effects break SUv(4?i/) x SUa(4?i/)) and too much (the target symmetry is SUv(n/) x SUa («/))■ 
Confusion stemming from too much symmetry can be avoided via the Gedanken algorithm dis- 
cussed in sect. ^ It explains why rooted theories have "extra" particles, in particular pseudo- 
Goldstone bosons, and why the extra symmetry protects physical (single-taste) correlators. Most 
of the criticisms of Refs. [ 12, |l^, 14, 15, |l^] can be refuted with this framework [ 17, 83]. 

Valid criticisms of rooted staggered fermions should focus on the difference between staggered 
fermions and the continuum, namely on the interactions that break the full SUv(4«/) x S\]\{Anf) 
taste-flavor symmetry. The most glaring issue is the violation of unitarity. With the full symmetry 
such violations cancel identically (for physical correlators). Without the full symmetry they do not 
cancel, but rooted staggered chiral perturbation theory (RS;^PT) offers a way to describe them. Fits 
to pseudoscalar masses and decay constants (light and heavy-light), and studies of the two-particle 
contribution to the scalar propagator, give evidence that RS;^PT works. It not only guides the chiral 
extrapolation but also provides a framework for estimating the associated uncertainty. 



Refutation of specific criticisms, here and in Refs. [17, 83], does not prove that rooted, stag- 
gered fermions are valid. It is remarkable, however, that so many numerical tests have shored up 
the theoretical framework, when any one of them could have gone wrong. Of special interest here 
is the interplay of chirality and topology. The analysis of chirality in Refs. [|12[ [l^, 15, 16] is 
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wrong, because it discusses a taste-nonsinglet chiral symmetry that has nothing to do with topol- 
ogy. The conventional staggered taste-singlet chirality agrees with the index theorem, in the same 
way any two lattice definitions of topological charge do. 

More tests can and should be carried out. A comprehensive study of the scaling of the Shamir 
defect, aA„, should help decide whether all taste breaking interactions vanish in the continuum 
Umit. For self-consistency, this study should be done with a rooted staggered sea. Tests of the 
sea, such as the scaling of the two-particle contributions to taste-singlet correlators, would bolster 
confidence in RS;tPT. Most of Creutz's discussion of 't Hooft vertices can be refuted with the setup 
in sect. Q but it does bring out the need to check numerically whether the members of near-zero 
mode quartets are all of different taste. 

Much of our structural, and not to mention practical, understanding of rooted, staggered 
fermions depends on RS;^PT. It is based on plausible arguments []64|], with one reservation. An 
original justification for ;^PT is that it gives the most general description of particle interactions, 
consistent with unitarity [93]. What happens when unitarity is lost? Can a non-unitary ;^PT de- 
scribe a non-unitary gauge theory? Is, for example, cluster decomposition enough [0]? This is a 
basic concern wherever partially quenched ;^PT is used. With other formulations of lattice fermions 
one could, if necessary, avoid it. With staggered fermions it seems essential. 

A friend of mine, who expects that Eq. (LI) is not valid, says that the burden of proof is on 
the staggered community. He is correct, of course, but only half correct. A "proof" is unlikely to 
be mathematically rigorous. Instead, as in all of lattice gauge theory, methods will be validated, or 
not, with a combination of theoretical framework and numerical simulation. It does not make much 
sense for proponents to prove to themselves that their methods are acceptable. Skeptics need to be 
engaged, examine the theoretical and numerical evidence in favor of rooted staggered fermions, 
think about the issues clearly, and state where the shortcomings lie. In other words, their job now 
is to be like Mike. Think it through and write it up! 
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